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I. INTRODUCTION 

The rare decay B —)• K *has recently been attracting much attention from both the 
experimental nun and theoretical [TOl - l^ sides due to the various observables associated 
with this decay that are susceptible to reveal new physics (NP). In particular, reference 
m has brought to light an overall tension between the Standard Model predictions and 
the experimental data and has suggested that a modihcation to the Wilson coefficients 
could resolve this tension. 

To investigate signals of NP, one usually focuses on the region of the spectrum away 
from and ijj' resonances where short-distance (SD) interactions, as represented by Figs. 


onances are subtracted from the dileptonic spectrum. However, a careful analysis of 
B —)■ observables should consider the long-distance effects of the resonances in 

the SD dominated region. In this paper, we take into account the narrow resonance effects 
in the nonresonance region when calculating the differential decay rate and the forward- 
backward asymmetry in this decay. In doing so, we use a Breit-Wigner model for the 
resonances with momentum-dependent decay constants [22]. We note that the latter model 
hts the data on photoproduction and leptonic width of -0 and -0' simultaneously [23], is used 
for exclusive B —)■ for the first time. The effects of broad resonances, using quark 

hadron duality, are considered in Ref. HEIEU. 

In a previous paper [2S|, we have computed the full set of 7 independent B —)■ K* 
transition form factors. At low-to-intermediate q^, we used light-cone Sum Rules with 
AdS/QCD distribution amplitudes fDAs) [26] . These DAs are derived from the holographic 
AdS/QCD light-front wavefunction for iF*[27l [28] . We have fitted these with form factor 
predictions at high q‘^ from lattice QCD . In this work, the same method for the derivation 
of the form factors with updated inputs (B meson decay constant fs and b-quark mass rrib) 
are used to calculate the differential branching ratio and forward-backward asymmetry in 
the decay B —)■ . 

We find that including the resonance effects improves the agreement of our predictions 
with the LHCb data [1] and the latest CMS data [9] on the differential branching ratio. 
As for ApB, it seems that the inclusion of the resonances hardly changes our prediction for 
dimuon mass squared below the hrst resonance. Finally, we find that a negative shift in the 


1(a) and 1(b), are dominant. Experimentally, the region around the above two res- 
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B K* 

(a) SD penguin diagram 



B K* 

(b) SD box diagram 


B K* 

(c) LD resonance diagram 


FIG. 1 Feynman diagrams of the principal contribntions to the B —)■ K* decay. 


Wilson coefficient Cg enhances the agreement with the data for the differential branching 
ratio and the Afb at q‘^ below the hrst cc resonance. 


II. DIFFERENTIAL BRANCHING RATIO WITH RESONANCES 


In onr previous paper |25], we calculated the differential branching ratio for B —)■ without 

considering the effects of resonances. The inclusion of V’ and 'ifj' resonances, as illustrated 


in Fig. 1(c), is obtained by modifying Cg® with an additional term Gg®® which, using the 
Breit-Wigner model, can be written as 


Gr = (3Gi(/r) + G2(/x)) 
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( 1 ) 


where Gi and C 2 are the Wilson coefficients corresponding to the current-current operators 
Oi and O 2 evaluated at scale ^ mb and and F^(') are the decay constant and total 
width of the cc resonance respectively. We use the same convention for the effective 
operators as in reference [31] and the following dehnition for vector meson decay constant: 


{0\cj^c\V) = fvCf, . (2) 

Since ip and ip' resonances are off mass-shell for different from in i? —)■ K *, 
we need to consider the g^-dependence of their decay constants [22] : 


Mq^) = fv{0) 


{V = iP,iP') 


(3) 
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V 

fv(0) 

/y("iy) 

Cy 

dy 


0.54 

1.25 

0.54 

0.77 


0.043 

1.04 

0.043 

0.043 


TABLE I Parameters (in GeV-based units) used in the q‘^ evolution of fv- 


with the h function being related to the imaginary part of the quark-loop diagram: 


h{q^) 


1 

Ibvr^r 



20r 

IT 


+ 4(l + 2r) 



1 

- arctan 

2 



(4) 


where r = q^/Am?^ for 0 < rUg is the effective quark mass and assuming that the 

vector mesons are weakly bound systems of a quark and an antiquark, we take = my/2. 
As a result, Eq. dehned for 0 < < dm^, is an interpolation of fy from the experimental 

data on /y(0) (from photoproduction) and /y(my) (from leptonic width) based on a quark- 
loop diagram. We assume /y(g^) = fvi'^v) 4"^® numerical values of the 

parameters cy and dy in Eq. are given in Tab. Hira. 

The resonance contributions Gg®® augments the short-distance contributions in the 

effective Hamiltonian: 


= Cf - . (5) 

The minus sign in Eq. [^is due to our choice of convention for the Wilson coefficients. The 
real and imaginary components of Gg°* as a function of g^ are shown in Fig. To calculate 
the differential branching ratio including the resonance contributions, one should replace 
G|® by Gg®* in the differential branching ratio expression given in Ref. |25j . 

As for the 7 form factors which parametrize the B —)• K* transition, they are calculated 
using AdS/QCD DAs [2S] in conjunction with light-cone sum rules at low to intermediate 
g^. For high g^ values, we use the latest lattice data for B —)■ K* transition form factors 
[32j. Note that we use the lattice results reported under ensemble f0062 as they correspond 
to hner lattice spacing. We use the following two-parameter form to ht the form factors 
obtained from AdS/QCD at low-to-intermediate g^ and the lattice data at high g^: 


Hi) = 


F(0) 


+ 6 ; 


( 6 ) 


The updated values for G(0), a and b are given in Tab. Our prediction for the differential 
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FIG. 2 Plots of 3fJ(G9°*) and Q^(G 9 °*) versus q^. In the left plot, the solid curve is 
while the dotted curve is In the middle and right hgures, the solid and 

dashed curves correspond to utilizing momentum-dependent and momentum-independent 

decay constants, respectively. 


F 

Ao 

Ai 

A 2 

Ti 

T 2 

T 3 

V 

m 

0.243 

0.244 

0.244 

0.258 

0.239 

0.157 

0.297 

a 

1.618 

0.586 

1.910 

1.910 

0.525 

1.147 

1.934 

b 

0.561 

-0.356 

1.498 

1.082 

-0.459 

-0.114 

1.089 


TABLE II Updated £t parameters for the seven independent B —)■ K* form factors used in 

Eq. g 


branching ratio including the effects of the resonances ip and ip' as obtained by using the 
above form factors is shown in Fig. where we compare with the latest data from LHCb 
[1] and CMS [9]. Our numerical results are calculated with the input parameters given in 
Tab. |TTT] and the Wilson coefficients tabulated in Tab. W Fig. clearly shows the effect 
of including resonances with the momentum-dependent decay constant on our prediction of 
the differential branching ratio. 
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FIG. 3 The AdS/QCD prediction for the differential branching ratio of the 
B —)■ K*iji~ decay, with (solid red) and without (dashed blue) resonances, as compared 
with the latest LHCb —)■ (diamonds)) and CMS (5° —)■ (crosses)) 

data. Note that this plot is qualitative and our predictions for each experimental bin for 
this observable are shown in Tab. |V| in Appendix A. 


III. FORWARD BACKWARD ASYMMETRY 


The forward-backward asymmetry distribution in dileptonic rare B —)■ decay is 

defined as: 


dApB _ 1 ( [^ . . d^T 

dq^ dT/dq'^ yJo ^ dq^d cos 61 


/ Ofd) —--- 

/_! dq^dcos&e 


( 7 ) 
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where Qg is the angle between the positive muon and the line of flight of K* in the p"''/! 
rest frame. This distribution to next-to-leading order(NLO) accuracy in is given by 


dA 


FB 


GfirrVttI 


dg2 





1287r3 

^tot ^ 

47r 


Op 


dvr 




, , , ,T2(g2^ 

[vriB rriK*) .. + (ms - m^.)- 


{niB + rriK*) 


V{q^) 

1 

vW) 

-1 




X 


Cf + 


dvr 


-c'r Hr) 


(ms - niK*) 1 - ^ 


m 


B 


7li(g2 


X 


agCp TT^ fBfK*,±^B,+ 


dvr iVp 


dtid) 


ms 




»T 


(nf) 


M 


( 8 ) 


where $k*,± is the transverse twist-2 DA for K*. The NLO contribution in Eq. |^is directly 
sensitive to this DA and therefore it would be interesting to examine its relative signihcance. 

Our prediction for Ape distribution is given in Fig. in which the latest data points 
from LHCb, including the zero-crossing point Qq = 3.7l5i GeV^[8], and CMS [9] are shown 
as well. 


IV. RESULTS 

The AdS/QCD predictions for the B —)■ A'*/r+/i“ differential branching ratio are shown 
in Fig. 1^ We can see that the resonance effects are significant and improve the agreement 
with the experimental data for regions above m^. The gray bands in this hgure (and in all 
subsequent hgures) represent the uncertainty due to the renormalization scale p (taken in the 
range mb/2 < /i < 2mb) and the error bars on the lattice data for the form factors. The latter 
is dominated by the uncertainty in A 2 lattice calculations. Fig. [^a) shows our prediction for 
the differential branching ratio when we assume a momentum-independent decay constant for 
■0 and "0' (dashed curve). We note from this graph that the only signihcant difference occurs 
at below the first resonance. As is the case for the inclusive B —)■ assuming 

momentum-dependent decay constants leads to better agreement with the experimental data 
for small q^. Fig. [^c), on the other hand, shows our predictions for the differential branching 
ratio when additional NP contributions are added to the Wilson coefficients Cf^ and 
We note that assuming = —1.0 and = —0.01, as suggested by the authors of Ref. 















FIG. 4 LO (dashed) and NLO (solid) predictions for ApB inclnding (red) and exclnding 
(bine) the resonance effects. We compare to the latest LHCb (diamonds) and CMS 
(crosses) data. Note that this plot is qualitative and our predictions for each experimental 
bin for this observable are shown in Tab. 


VI 


in Appendix A. 


[33] . produces better agreement with the data, especially at high q^. In Fig.j^e), we compare 
our predictions with those obtained from sum rules (SR) DAs. It seems that AdS/QCD DAs 
produce results generally lower than those obtained from SR DAs [34] . especially for larger 
. The predictions for each experimental bin for this observable are shown in Tab. |^in the 
Appendix. 

Our predictions for Aps are shown in Fig. First, we observe that the leading-order 
predictions miss all but one of the experimental data points as well as the zero-crossing 
point. Second, as pointed out in Ref. |3T], the inclusion of NLO contributions leads to a 
signihcant shift to the zero-crossing point (of order 30%) and an overall better agreement 
with the most recent data on App below the hrst resonance. We observe that the inclusion 
of the two resonances does not have any noticeable effects for this observable outside the 
resonance regions. Consequently, as shown in Fig. [^b), assuming momentum-independent 
decay constants for -0 and ip' does not change our predictions signihcantly. On the other 
hand, assuming the NP contributions = —1.0 and = —0.01, produces much better 
agreement with the experimental data, as seen in Fig. |^d). Finally, predictions for App 
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q" [GeV"] 


(a) The differential branching ratio using 
(/^-dependent (solid) and q^-independant 
(dashed) fv- 



(b) Afb using g^-dependent (solid) and 
-independant (dashed) fv- 



[GeV"] 


(c) The differential branching ratio within 
the SM (solid) and with new physics 
(Cl^P.C^P) = (-1.0,-0.01) (dashed). 
The bin by bin predictions are given in the 
Appendix A. 



(d) Afb within the SM (solid) and with 
new physics = (-1.0,-0.01) 

(dashed). The bin by bin predictions are 
given in the Appendix A. 



q^ [GeV"] 


(e) The differential branching ratio using (f) Afb using AdS/QCD (solid) and SR 
AdS/QCD (solid) and SR (dashed) DAs. (dashed) DAs. 

FIG. 5 Variations in the AdS/QCD predictions of the differential branching ratio and Aps 
as explained in each figure caption. The red and blue curves show the results with and 
without the inclusion of if and if' resonances. 
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based on AdS/QCD DAs are more or less in similar agreement with the data as those 
obtained from SR DAs, as illustrated in Fig. |^f). The predictions for each experimental 
bin for this observable are shown in Tab. |VI| in the Appendix. 


V. CONCLUSION 

We used the form factors and DAs as predicted by the AdS/QCD correspondence as 
well as taken into account the possible cc resonance contributions to give predictions for the 
B —)■ K*ii^ differential branching ratio and forward-backward asymmetry. The inclusion 
of 'ip and Ip' resonances is done by using the Breit-Wigner model with momentum-dependent 
decay constants. This leads to better agreement with the experiment data for the differential 
decay rate outside the resonance regions. However, the forward-backward asymmetry outside 
the resonance region is not affected by the presence of resonances. We conhrm that a 
negative contribution to Cg and a small contribution to Cy, as suggested in Refs mm, 
leads to better agreement with the experimental data. Comparison of predictions from 
AdS/QCD DAs and SR DAs shows that the former produces better or identical results when 
compared with experimental data on the branching ratio and Aps- It would be interesting 
to investigate the use of our AdS/QCD form factors and DAs to compute other angular 
observables associated with B —)■ decay for the whole range of in particular, the 

observable for which there is a discrepancy between the theory predictions and the LHCb 
measurement [35] . 
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«s(mz) = 0.1185 

aem = 1/133 

mz = 91.19 GeV 

=0.119 

/x* = 0.225 

fs = 0.18 

MB(Borel) = 8 GeV 

So = 36 GeV 


niq = 0.35 GeV 

niB = 5.28 GeV 

rus = 0.48 GeV 

rriK* = 0.89 GeV 

me = 1.4 GeV 

mp = 3.10 GeV 

rrib = 4.6 GeV 

m^> = 3.69 GeV 

mt = 173.5 GeV 



TABLE III Numerical values of the input parameters used in our calculations. 


Cl 

C2 

Cs 

C4 

Cs 

Ce 



Cg 

Cio 

-0.148 

1.060 

0.012 

-0.035 

0.010 

-0.039 

-0.307 

-0.169 

4.238 

-4.641 


TABLE IV Values of the Wilson coefficients at /i = rrib. 


Appendix A: Numerical inputs and bin by bin results 


Throughout our analysis, we have used the input parameters presented in Tab. |III| where 
all quark, meson and the intermediate boson masses, as well as the experimental value 
of as{mz), are taken from the latest Review of Particle Physics [36]. The two K* decay 
constants, fx* and /^*, are AdS/QCD predictions which are dependent on the masses of 
the quarks in the K* meson |25j . 

We use the next-to-next-to leading order evolution for the strong coupling constant 
which can be found in the Appendix of Ref. EH- We also present the values of the 10 Wilson 


Coefficients at scale /i = ruf, in Tab. |IV| The complete set of equations used to obtain these 
values have been collected in the appendix of Ref. 
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Bin (GeV^) 


(dB^^^/dq^'^ 

(dB^P/dq^'^ 


Experiment 

Process 

[0.10-2.00] 

"J-' '^+0.076 

n 77S“0-044 
'O+0.075 

0 7(i7“0.050 

O'-‘ll'+0.078 

0.768;°:“^ 

0.5921°;}?? 

to 

+ 

i 

-X- 

+ 

"R 

+ 

•R 

1 

"P 

ffi 

0 

[2.00 - 4.00] 

n c;9Q+0.044 
U.JZO_Q 037 

n c; 97+0-045 
U.oz/_oo37 

U.+OZ_o 045 

0.4551^^? 

0.5591°;}“ 

[4.00 - 6.00] 

U.DD^_q Qg3 

U.DDO_o 063 

U.OOZ_o 068 

0 444+0-072 

U.OOO_o 080 

0.2491°;}?? 

[6.00 - 8.00] 

0.869l[];i°? 

u.aau_o 105 

u./ua_0 096 

D 744+0+12 

U. 1OO-0 092 

f) 440+0-136 

U-OOU_0 123 

[11.00- 12.50] 

1.2861°;^^' 

1 174+0-128 
1-1 '+-0.121 

1 044+0-137 
r-U+O_0 130 

f) 966'*'°'^^^ 
U.yDD_o 081 

n S9S"I"9-214 
U.O^O 0.i97 

[15.00 - 17.00] 

i.OfU_o 123 

1 108+0+16 
t-ty8 o 105 

1 114+0 145 
1-11+-0.138 

0.990l°:Ji? 

n 644+°-“3 

U.D++-0 159 

[17.00 - 19.00] 

i.Uf Z_o 125 

0 084+0-029 
U.a54_0.ii4 

n C79+0-027 
U.5/2 ;o 050 

0.8071°:°?? 

f) 110+0-099 

U.iiD-o 084 

[1.00 - 2.00] 

U.OiU_|_o 010 

0 51 2“0-005 
u.diz_|_o 010 

n 47Q-0.009 

0.^:10^0.024 

n 474-0.008 
'-’-+‘'3+0.023 

0 47+0.076 
'-’-+‘-0.076 

ta 

0 

i 

* 

0 

■R 

+ 

■R 

1 

0 

L5 

[2.00 - 4.30] 

n c;q 9+0.046 
U.JOZ_o 040 

0 >54I5+0-047 

U.090_o 040 

n ac;q+o.o 2 o 
U.455_Q Q 3Q 

0 461+0-020 
U.+Di-o 030 

n qq+O-045 
'^•'^^-0.045 

[4.30 - 6.00] 

0 f)7fi+0-065 
U.Dm_o 053 

0 f)77+0-065 
0'-'-’' ‘ -0.053 

0 561 +0-052 

U.ODi_o 045 

U.ODZ_o 045 

n Q/1 +0-058 
'^•'^^-0.058 

[6.00 - 8.68] 

u.ai^_o 100 

1 qoi+0.120 
l-U9i_o 115 

n 744+0.089 

O’-‘+'3-0.082 

n cq^+O-095 

U.094_o 082 

n A '7+0.050 
'^•^^-0.050 

[10.09 - 12.86] 

r-oau_o 121 

1-148^°;^°^ 

1 no7“l“0.106 
i.UZf _o 091 

f) 922+0-078 
u-y^^0.063 

n (qo+0.064 

'^•'^"^-0.064 

[14.18-16.00] 

1 1f)4+0-300 

i.iD4_0 115 

1-0581°:?“ 

1 1 qc;+0.027 
'^^-0.089 

o.96o1°:°?9^ 

f) 05+0-078 
'-’-'3'3-0.078 

[16.00 - 19.00] 

1 1f)4+0-300 
i.lD4_0 115 

1.0581°;?“ 

n Q4fi+0 027 

U.y4:D_o 089 

0.8681°:°?? 

0 42+0-042 
'3-+^-0.042 


TABLE V Bin by bin valnes of the branching ratio dehned as x /g? ^i] 

bin, with and withont resonances as well as new physics contribntions as compared with 

the latest LHCbP and CMS [9j data 
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Bin (GeV^) 

(^Afb^ 

{^Tb) 

{a%) 


Experiment 

Process 

[0.10 - 0.98] 

u.uaD_,_g Q^Q 

-0.095;[];™3 

-0.103;°:0?5 

-0.103;°:“^ 

—0 003+'^-''®® 

U.UUO_o 060 

ba 

o 

* 

o 

7: 

+ 

■R 

1 

'p 

o 

[1.10 - 2.50] 

u.uya_g Q^g 

-0.097l[];°li 

—0 140+*^-°^° 
U.l+U_o 020 

n 1 QQ+o.oio 

U.iO5 0 021 

-0 iQi+O-070 
u.ryi-o 079 

[2.50 - 4.00] 

U.Uii_o 005 

U.UiU_o 005 

U.UDO_o 017 

U.UD.5_0 017 

n 1 1 O+0.075 

U.118_o 088 

[4.00 - 6.00] 

u.ur0_o 020 

0 n7i5+0-005 

u.uro_o 020 

U.UZO_o 009 

U.UZO_o 007 

U.UZO_o 050 

[6.00 - 8.00] 

029 

U.i+D_,_o 031 

0 104+0009 

U.iUO_o 018 

nil O+0-006 
U.iiO_o_oi9 

U.lOZ-0 041 

[11.00 - 12.50] 

+0.012 

n 1 Q9~0.015 

U.lOZ_^0.020 

f) 1 c:o+0.005 
U.iOO_o 005 

U.i+U_o 005 

n 91 O+0-041 

U.ol5 o o41 

[15.00 - 17.00] 

n 1 o-| —0.015 
U.lol^0.019 

U.iDO_,_o 017 

0.104;[];E]J2 

n 1 9 o~0.008 
U.1ZO_^0.009 

n 411+0-041 
'-'•+11-0.036 

[17.00 - 19.00] 

U.1ZO_,_0 012 

+0.011 

0.082;[];E]io 

u.uyi+o 008 

f) 4n5'’-0-049 

U.OUO_o 050 

[1.00 - 2.00] 

U.ii+,o 016 

U.liO+o 016 

_n it;n-0 007 

U.10U_,_o 010 

n 1 T q~0-009 
11.145^0.011 

n 97+0.184 
-0.406 

ba 

o 

i 

* 

o 

7: 

+ 

■R 

1 

+ 

L2 

[2.00 - 4.30] 

-o.oi8iO:;]i° 

-0.017l[];[]J0 

—0 071+0 013 
'-'■O’' 1-0.017 

—f) 069^®'^^^ 

u.uoy_o 016 

n 1 9+0.158 
'^‘-'■^-0.149 

[4.3 - 6.00] 

0-08ll°;“^ 

0.08ll‘]:°°^ 

U.UOZ_o oil 

u.uaz_o on 

f) 03+0-153 
O’-'-' '1-0.153 

[6.00 - 8.68] 

U.i+y_o 029 

U.iO+_o 030 

nil q+O-006 

U.ilO_o 022 

U.lzu_o 027 

0’-'-’+-0.101 

[10.09 - 12.86] 

n 1 Qri~o.oii 

u.iyo,o 012 

U.iOi+o 005 

0 103+0011 
U.iDO_o 014 

0 lOQ+O-Ooe 
u.iuy_o 019 

f) IR+0-061 

0’-l'l-0.061 

[14.18 - 16.00] 

u.iy ^_,_0 019 

U.lOO+o 014 

0 IfiQ+O.OOS 

U.iDO_o 010 

0 lOQ+O-OOS 

u.iuy_o 007 

0 40+0-041 
0’-+0’-0.061 

[16.00 - 19.00] 

0 i4n“0.oi3 
U.i+U_,_o 014 

U.lOi+o 013 

nil Q—0.010 

U.118^0.010 

U.lUZ+o 007 

n 33+0-071 
I’-'l'1-0.071 


TABLE VI Bin by bin values of Aps defined x J ^2 dq^i with and without 

resonances as well as new physics contributions as compared with the latest LHCb j8] and 

CMS [9] data. 
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